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Abstract 

Wc reexamine the constrained version of the Next-to-Minimal Supersymmetric Standard 
Model with semi universal parameters at the GUT scale (CNMSSM). We include constraints 
from collider searches for Higgs and SUSY particles, upper bound on the relic density of dark 
matter, measurements of the muon anomalous magnetic moment and of B-physics observables 
as well as direct searches for dark matter. We then study the prospects for direct detection of 
dark matter in large scale detectors and comment on the prospects for discovery of heavy Higgs 
states at the LHC. 

1 Introduction 

The Next-to-Minimal Supersymmetric Standard Model (NMSSM) is a simple extension of the 
MSSM that retains two of its most attractive features - a solution to the hierarchy problem and 
a natural dark matter (DM) candidate - while solving the naturalness problem. This is achieved 
by the introduction of a gauge singlet superfield, S. The vev of this singlet field determines the 
effective /x parameter which is then naturally of the EW scale [1-10]. In the NMSSM the upper 
bound on the lightest scalar Higgs mass can be larger than in the MSSM [11-15] or below the 
LEP bound if the latter decays into two lighter pseudoscalar Higgs states [15-20]. In addition, a 
very light pseudoscalar Higgs state is not excluded by LEP or B-physics observables [21-24]. The 
neutralino sector of the NMSSM contains an extra state, the singlino, which couples only weakly to 
most of the non-singlet (MSSM like) particles. Owing both to modifications in the Higgs sector and 
the neutralino sector, the DM properties can differ from those of the MSSM [25- 31] . In particular a 
LSP with a large singlino component has different annihilation properties than the bino LSP that is 
in general found in the CMSSM. It is even possible to have a very light singlino LSP if accompanied 
by a very light scalar. Even when the LSP has no singlino component, the more elaborate Higgs 
sector of the model provides additional channels for rapid annihilation through Higgs exchange. 
This can have implications for direct and indirect detection rates [31-35]. 

When imposing GUT scale boundary conditions, the constrained NMSSM (CNMSSM), contains 
only a handful of parameters, making it convenient for phenomenological analyses. In our definition 
of the CNMSSM, we impose universal parameters at the CUT scale except in the singlet sector. 
The model thus contains 6 free parameters plus an arbitrary sign [36]. Because of the relations 
among the various physical parameters at the weak scale, there are stringent constraints on the 
fundamental parameters of the CNMSSM at the GUT scale. Constraints from Higgs and SUSY 
searches at colliders as well as cosmological constraints, resulting from the precise measurement of 
the relic density of DM were examined in ref . [30] . Regions of parameter space where the predictions 
for DM differed from the ones of the CMSSM were highlighted. Apart from the additional allowed 
parameter space due to relaxed constraints from LEP on the Higgs sector, it was found that rapid 
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annihilation regions (where the LSP annihilates efficiently through a Higgs resonance) could occur 
at small values of tan/3, a direct consequence of the more elaborate Higgs sector of the model. 
However the scenario with a very light singlino LSP annihilating through a very light (and mainly 
singlet) scalar Higgs state, that can occur in the general NMSSM [28, 29], cannot be realised in the 
CNMSSM. 

Within the hypothesis of minimal flavour violation, the B-physics observables put additional 
constraints on the parameter space of the model [21-24,37], in particular at large tan/3 and when 
the heavy Higgs states are not too heavy. Explaining the discrepancy with the SM expectations for 
the muon anomalous magnetic moment also restricts the parameter space, in particular sfermions 
must be light [24, 38, 39]. Finally the experiments that have been searching for elastic scattering of 
DM in large detectors have now reached a sensitivity level where they can constrain the parameter 
space of SUSY models including that of the (C)MSSM [40,41] and the (C)NMSSM [31,33,34]. In 
particular the LSP with a large higgsino component that is found at large niQ and large tan/3 can 
lead to cross sections above the present limit. This can also occur in models with a rather light 
heavy Higgs doublet. 

Our first goal in this paper is to reexamine the parameter space of the CNMSSM including 
in addition to the cosmological constraints already considered in ref. [30], the constraints from 
B-physics observables [23,42,43], from the muon anomalous magnetic moment [38,44] and from 
the recent upper limit on DM direct detection [40,45]. After having identified the allowed regions 
we concentrate on those which differ markedly from the MSSM, where either the LSP has some 
singlino component or there is a possibility of LSP annihilation through a (singlet) Higgs resonance. 
We then explore the potential of future direct detection experiments for probing the remaining 
parameter space of the model. We find that with the exception of the singlino LSP scenarios, large 
scale detectors will probe most of the parameter space of the model. Signals are in several cases 
expected just beyond the present limits. 

As a secondary goal of this paper we address the issue of distinguishing the CNMSSM from the 
CMSSM in regions favoured by the cosmological measurements. Assuming that a signal compatible 
with SUSY is observed in flavour physics, DM detection and/or at the LHC, we investigate the LHC 
potential for probing the NMSSM Higgs sector. The Higgs sector of the NMSSM includes, as in 
the MSSM, a light scalar state h, two almost degenerate heavy scalar /pseudoscalar states H/A, 
and acharged state H^. In addition, the NMSSM Higgs sector contains two extra singlet states, 
a scalar ,s and a pseudoscalar a, whose mass can be below or above that of the (MSSM like) 
doublet states. A clear indication in favour of the NMSSM would therefore be the discovery of at 
least three neutral states with large mass splittings. The difficulty in observing the extra singlet 
states, s or a, is a consequence of their tiny couplings to SM particles unless their mixing with the 
doublet states are substantial. The LHC potential for probing the Higgs sector of singlet extensions 
of the MSSM, including the NMSSM, was explored in [17-19,46-51] with emphasis on the light 
Higgs doublet/singlet. The production of the heavy states was not considered. The possibility 
of distinguishing the NMSSM from the MSSM by observing a singlino LSP will not be discussed 
here [52-57]. 

The paper is organized as follows: in sec. 2 wc briefly describe the model. In sec. 3 we present the 
various constraints on the model. Typical case studies are presented in sec. 4 The direct detection 
prospects are discussed in sec. 5 followed by Higgs searches at the LHC in sec. 6. Our results are 
summarized in sec. 7. 

2 The CNMSSM 

We discuss here the NMSSM with the scale invariant superpotential 

W = XSH^Hd +'^S^ + (Yukawa couplings) , (2.1) 
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where no supersymmetric dimensionful parameters as ^ are present in the superpotential, i.e. the 
weak scale originates from the soft SUSY breaking scale only. The soft SUSY breaking terms in the 
Higgs sector are then given by 

Koft = rn%jHu\^ + mjjjH£ + ml\S\^ + (^XAxHuHaS + ^kA^S^ + h.c.^ . (2.2) 

A possible cosmological domain wall problem caused by a global symmetry [58] can be avoided 

by introducing suitable non-renormalizable operators [59, 60] which neither generate dangerously 
large singlet tadpole diagrams [3, 61-63] nor affect the low energy phenomenology. Other models 
solving the MSSM /x problem with extra gauge singlet fields and with no domain walls include 
the nearly minimal SUSY model (nMSSM), with an additional tadpole term for the singlet in the 
superpotential and/or in the soft scalar potential [64-66], and the UMSSM, with an extra U^l)' 
gauge symmetry [15, 67-71]. The DM properties as well as the possibility of generating the baryon 
asymmetry of the universe within these models have been analysed in refs. [34, 72-81]. 

We consider here a constrained version of the NMSSM (CNMSSM) with semi-universal param- 
eters defined at the GUT scale. Embedding the (N)MSSM in an underlying theory and assuming 
some spontaneous symmetry breaking mechanism (for example supergravity) typically leads to 
simple patterns for soft SUSY breaking parameters at the GUT scale, hence reducing significantly 
the number of arbitrary parameters in the model. The choice of semi-universality is motivated by 
the fact than imposing strict universality makes it difficult to find a consistent set of parameters 
in the NMSSM [82]. The reason is the following: in the CMSSM, in addition to the soft terms 
AqjUIq, Ml/2 defined at the GUT scale, Mz and tan/3 at the weak scale are used as inputs. The 
two minimization equations of the Higgs potential w.r.t. the two real Higgs vevs hu and ha are 
used to compute fi and B in terms of the other parameters (this leaves the sign of /i as a free 
parameter). Both fi and B have only a small effect on the renormalization group equations (RGEs) 
of the other parameters (in numerical codes, this is usually solved by an iterative procedure). On 
the other hand, in the CNMSSM neither /x nor B are present and one has to cope with three 
coupled minimization equations w.r.t. the Higgs vevs h^, and s. This means that starting from 
strict universality, with Mi m,Q, Aq as well as A and k as free parameters at the GUT scale, one 
usually ends up with the wrong value of Mz at the weak scale (as no dimensionful parameter is 
left to tune the correct value) . In addition, tan/3 cannot be a free parameter in this approach. 

The semi-universal approach described here was first used in the program NMSPEC [36]: First, 
let us assume that A as well as all the soft terms (except m|) are known at the weak scale (e.g. 
after integration of the RGEs down from the GUT scale,). One can define effective (s dependent) 
parameters at the weak scale: 

/X = As , v = Ks , B = A\ + i'. (2.3) 

The minimization equations w.r.t. hu and hd can then be solved for the effective /x and B, as in the 
MSSM, in terms of the other parameters (inch Mz and tan/3). From jj, and B one then deduces 
(for A and A\ given) both s and k. Finally, from the minimization equation w.r.t. s, one can easily 
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equations giving /x (up to a sign), B and m| read: 



obtain the soft singlet mass mi in terms of all other parameters. At tree level, the minimization 



,,2 ^^d ffu _ l\^Z 



mjf^ — mjf^ tan^/3 
tan2/3 - 1 



B = '^{m%^+m%^ + 2^^^ + X\hl + hi)), (2.4) 

m| = X^^{A, + 2i.)-,yiA, + 2,y)-X\hl + hl). 
A* 

The radiative corrections to the scalar potential show a weak dependence on s, k and m| which 
can be included in the minimization equations. These become non-linear in the parameters to 
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ml = { 008^2/3 + ^ sin22/3 ) M| , (2.5) 



solve for and have therefore to be solved iteratively. The derived parameters k and nig affect the 
RGEs of the other parameters not only through threshold elfects around Mgusy; but also through 
the (3 functions. However, the numerical impact is relatively small such that an iterative procedure 
converges quite rapidly again. 

As rrvg is computed from the minimization equations, it is difficult to find parameters such that 
it assumes the same value as the Higgs doublet (or other scalar) soft masses squared at the GUT 
scale. On the other hand, the mechanism for the generation of soft SUSY breaking terms could 
easily treat the singlet differently from the other non-singlet matter multiplets [83]. This would 
also affect the coupling which involves the singlet only. Hence, in our semi-universal approach, 
K, s and m'g are computed from the minimization equations leaving the following free parameters: 

• tan/3, sign(//) at the weak scale; 

• A at the SUSY scale; 

• Ml/2, 1^0, Aq and A,^ at the GUT scale. 

Not all choices of parameters are allowed in the CNMSSM. Some lead to negative mass squared for 
scalar fields (Higgs or sfermions), others to Landau poles below the GUT scale for the dimensionless 
couplings. The Landau pole condition leads to an upper bound on A<.7 which depends on tan/3. 
The LEP bounds on the Higgs mass further constrain values of A> .L This is due to the the fact 
that although the diagonal entry for the light double h in the scalar Higgs mass matrix increases 
with A: 

,2 _ / „„„2o« , ^ „;„2o« \ 71^2 

g' 

the mixing between the doublet and the singlet states is proportional to A times some combination 
of soft SUSY breaking terms. In the general NMSSM, this combination can be made to vanish. 
In the semi-universal NMSSM however, this is no more the case. Therefore, if A is too large 
(A > .1) mixing effects lead to a light eigenstate in the Higgs spectrum with a substantial doublet 
component. Such a state is usually excluded by LEP, unless it decays to two lighter pseudoscalar 
singlet Higgs states. Such a situation requires however some fine tuning, especially on A^, and will 
not be studied here. 

When A is small the mixings between the Higgs doublet and singlet states are small. As in the 
MSSM, the masses of the nearly degenerate heavy Higgs doublet states read 

2 2 

We will denote the scalar singlet s and the pseudoscalar singlet a and their masses read, respectively: 

ml = v{A^ + Au) , ml = -ZvA^ . (2.7) 

The parameter being only slightly renormalized from the GUT scale down to the SUSY scale, 
eq. (2.7) shows that the masses of the singlet states are proportional to the value of A^ at the GUT 
scale. The condition that both squared masses are positive together with eq. (2.3) implies 

-A{B-Axf < A^{B-Ax) < 0. (2.8) 

The parameter B, obtained from the minimization equations (2.4) depends on mo, M1/2, A and 
tan/3, while Ax depends on Ao, A and tan/3 through the RGEs. This means that, for sign(;u) positive 
(which we will always assume in the following), either A^. > and Aq > Ao{mo, M1/2, X,tan.p) or 
Afi < and Aq < Ao(mo, M1/2, A, tan/3). Moreover, for A^ > and tan/3 moderate, large values of 
mo or Ml/2 (implying B large and positive) lead to a negative mass squared for the pseudoscalar 
singlet a and are therefore disallowed. If Af^ > and tan/3 is large, however, B remains small up 
to large values for mo and M1/2 which are no longer excluded. 
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3 Constraints 



To evaluate numerically the supersymmetric spectrum we used the NMSPEC program from the 
NMSSMTools package [36]. Starting from GUT scale parameters and using the RGEs, this code com- 
putes the Higgs spectrum including higher order corrections [12, 13] as well as the masses of sparti- 
cles at one-loop. The NMSSMTools package also includes all the available experimental constraints 
from sparticle and Higgs searches at LEP and Tcvatron (for details on the exclusion channels see 
refs. [36, 84, 85]) as well as constraints from B-physics [23] and from the muon anomalous magnetic 
moment [38]. For the computation of the relic density of DM we rely on micrOMEGAs [86] which 
is included in the NMSSMTools package. The elastic scattering neutralino nucleon cross section is 
computed with micrOMEGAs2.2 [87]. We used the central value of the top quark mass measured at 
the Tevatron, mt = 172.6 GcV, and chose mb{mb) = 4.214 GeV and as{Mz) = .1172. Note that 
the top quark mass does affect the value of the light Higgs mass as well as the calculation of the 
spectrum, especially at large mo. As stated above, we assumed sign(/i) > 0. 

3.1 B-physics 

Assuming minimal flavour violation, a severe constraint on the CNMSSM originates from the 
branching ratio Br(b sj). We require that the theoretical prediction, including uncertainties, 
falls within the 2a range [43]: 

3.07 X 10-'' < Br{b s^) < 4.07 x 10"^ . (3.1) 

The most important SUSY contributions arises from chargino-squarks loops and charged Higgs- 
top quark loops. For > considered here, the chargino exchange diagram gives a negative 
contribution relative to the SM one so that the branching ratio for b — > 57 drops below the allowed 
range when the chargino is light, that is in the low mo, Mi^2 region. As in the CMSSM there is a 
strong dependence on At (hence ^0) from the mixing in the stop sector. The chargino and charged 
Higgs contribution partially cancel each other when At > thus the b ^ sj constraint is most 
effective for negative values of At (or Aq). Also as in the CMSSM, this constraint is stronger for 
large values of tan/3. 

The neutral Higgs contribution to the branching ratios for Bg and the charged Higgs 

contribution to -B+ — r'^Uj- both feature a strong tan/3-enhancement (respectively oc tan^/3 and 
tan^/3). The latest 95%C.L. bound from CDF at Tevatron [42], 

Br{Bs /i+Ai~) < 5.8 x 10"^ (3.2) 

and the 2a world average [43] 

0.34 X 10"^ < Br{B+ t+u^) < 2.30 x 10"^ (3.3) 

thus constrain the large tan/? regions where the heavy neutral/charged Higgs are rather light. These 
conditions arc met at low values of m,Q as well as in the so-called 'focus point' region at large mo- 
The relative importance of each channel is influenced by the trilinear coupling At which control 
the flavour violating neutral Higgs coupling. In particular for Bg — > chargino/stop 
contribution is suppressed for small values of At (obtained when Ao is negative). On the other 
hand S+ t'^Ut is independent of the trilinear coupling while negative values of ^0 strenghten 
the constraints from 6 ^ 57. 

The constraints from AMs,AM(i are also taken into account in NMSSMTools [23]. Given the 
present experimental results these observables do not allow to constrain additional parameter space 
relative to other B-physics observables. 
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3.2 Muon anomalous magnetic moment 

The anomalous magnetic moment of the muon shows a deviation from the SM prediction when 
using only e^e~ data to estimate the leading order hadronic SM contribution [44]. We imposed 
the 2(7 bound: 

8.7 X 10~^° < (5a^ < 4.6 x lO"'^ . (3.4) 

As in the CMSSM, the most important supersymmetric contribution comes from the chargino/sneu- 
trino exchange and is proportionnal to the Yukawa coupling (which means oc tan/3). To explain the 
deviation from the SM prediction this observable therefore requires light charginos and smuons, 
hence favours the low mo, region. At larger values of tan/3 the allowed band moves towards 

larger values of mo, since light sparticles can then give a contribution that is too large [38]. 

3.3 Relic density 

The allowed range for the DM relic density from cosmological data [88, 89] is 

.094 < D.h^ < .136 (3.5) 

when including some uncertainty on the cosmological model parameters [90]. We only use the 
upper bound on Clh^ since we assume that the neutralino does not necessarily account for all the 
DM. 

A complete calculation of the neutralino DM in the (C)NMSSM [28-31] has shown that often the 
same mechanisms as in the MSSM for neutralino annihilation are at work. Nevertheless important 
differences are found: the extra pseudoscalar singlet Higgs state a opens up the possibility of 
resonant annihilation and the singlino component of the neutralino LSP can alter the prediction 
for annihilation cross sections. The new mechanisms that can provide a DM candidate compatible 
with the WMAP results are: 

• Annihilation near the pseudoscalar singlet a resonance, occurring for A ~ .1. 

• Coannihilation of the singlino LSP with a stau or stop NLSP for 1 and small mo. 

• Coannihilation of the singlino LSP with a higgsino NLSP for A^ 1, large mo and tan/3. 

• Coannihilation of the singlino LSP with the bino NLSP, the bino rapidly annihilating 
through a Higgs resonance, for A^ 1 and large tan/3. 

3.4 Direct detection 

The best limits on the spin independent (SI) neutralino-proton cross section have been set recently 
by Xenon [45] and CDMS [40] . We will use the mass dependent CDMS limit since this detector has 
the best sensitivity for masses above ~ 50 GeV. In our model a lighter neutralino is excluded by 
the LEP limit on the chargino mass because of the universality condition on the gaugino masses. 
The CDMS limit corresponds to cTp^ > 4.6 x 10"* pb for m^ = 60 GeV and cTp^ > 2.5 x 10"'^ pb for 

= 670 GeV. Those limits are obtained by making standard assumptions about the DM density 
and velocity distribution as well as for specific choices of nuclear form factors. The theoretical 
uncertainties associated with the choice of the velocity distribution have been estimated to be about 
a factor 3 [91]. These uncertainties will not explicitly be taken into account here. When deriving 
constraints on the parameter space additional uncertainties that arise from our lack of precise 
knowledge of the quark coefHcents in the nucleon are taken into account by choosing two different 
sets of coefficients: the default values of micrOMEGAs [87] corresponding to {cr^j^, (Tq) = (55, 35) MeV 
and (45,40) MeV. The latter implies a lower s quark coefficient in the nucleon. First lattice 
calculations have recently obtained results that favour a not too large value for (TttAt [92]. 

The SI interaction is usually dominated by Higgs exchange provided the neutralino LSP has 
enough higgsino component to couple to a Higgs. The light Higgs doublet often gives the largest 
contribution although the contribution of the heavy Higgs doublet can dominate at large tan/3 when 
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its couplings to neutralinos and quarks are enhanced. As in the MSSM we therefore expect large 
cross sections when the heavy Higgs is rather light. This occurs both at low mo and in the focus 
point region at large mo. In the latter case the elastic cross section should be specially enhanced at 
large tan/3 due both to increased Higgs couplings and to a LSP with a large higgsino fraction. The 
squark exchange diagram also contributes but is suppressed by the large squark mass except in the 
low mo, Ml/2 region. This contribution is however never large enough by itself to obtain a^p near 
the experimental bound. Note that a pure singlino LSP has no coupling to the Higgs doublets nor 
to the squarks so one expects an extremely small SI cross section in this scenario. 

A higgsino LSP also couples to the Z boson and thus contributes to the spin dependent (SD) 
cross section. The experimental limits (Up^ ~ 10"^ pb [93,94]) are still orders of magnitude too 
large for this process to put constraints on the model. We will not consider it in the following. 

4 Results 

We first performed a general scan over the parameter space of the CNMSSM for fixed values of the 
SM parameters in order to find regions satisfying all theoretical and experimental constraints. We 
did not find any scenarios satisfying all constraints for small values of tan/3 < 5. This is mainly a 
result of the LEP constraints on the lightest Higgs boson coupled with the (g — 2)^ constraint. In 
particular, the allowed regions at small values of tan/3 = 2 analyzed in ref. [30] cannot explain the 
discrepancy in the measurement of the {g — 2)^ (this constraint was not included in this former 
analysis). We then picked fixed values for A and tan/3: A = .1 or .01, tan/3 = 5, 10 or 50 and for each 
choice scanned randomly on the remaining free parameters, namely mo, Mi/2i a-^id A^. As in our 
previous analysis, the regions in the parameter space allowed by all theoretical and experimental 
constraints correspond to either < and < Ao{X,tan/3) ot A^, > and > Ao(A,tan/3), 
the latter case appearing only if A = .01 and/or tan/3 = 50 (the other values of A, tan/3 always lead 
to light states in the Higgs sector excluded by LEP when A,^ > 0). 

Finally, we identified the values of Aq and A^, for which the main neutralino annihilation channel 
is the pseudoscalar singlet a resonance, or the LSP is mainly singlino. This will be illustrated in 
the following when we present plots in the mo,Mi/2 plane for selected values of A, tan/3, Aq and 
A,. 

4.1 Large A: singlet resonances 

When A is large ( > .1), the singlino component of the LSP is always small. Scenarios characteristic 
of the CNMSSM are therefore only those where the bino LSP annihilation proceeds through the 
pseudoscalar singlet a resonance. As explained in sec. 2, the larger A, the stronger the constraints 
from Higgs searches at LEP. These can still be evaded if is such that the light Higgs doublet 
h decays mainly into two light pseudoscalar singlets, h aa [19]. This requires however some 
fine-tuning on A^^. Furthermore, the a being very light, it cannot be used as a resonance in 
order to obtain the correct value for the DM relic density. In such case, one has to rely on the 
coannihilation of the bino LSP with light sfermions (occurring at small m,o) in order to satisfy the 
WMAP constraints, as in the CMSSM [19]. We will not discuss this case here. Another possibility 
is to maximize the Hght Higgs doublet mass by taking a small value of tan/3 (e.g. for tan/3 = 2 
and A = .5, as in ref. [30]). However, these regions are incompatible with the {g — 2)^ constraint 
as stated above. We therefore consider only A = .1 and tan/3 = 5, 10 and 50. In these cases, LEP 
constraints on the light Higgs doublet imply small < and large Aq < 0, except for tan/3 = 50 
where one can also have small > and large ^o > 0- 

As an example of a scenario with bino LSP annihilation through the a resonance we consider 
in fig. 1 the case ^o = -900 GeV, A^, = -60 GeV with (a) tan/3 = 5, (b) tan/3 = 10. The a 
resonance is found around M1/2 ^ 350 GeV, ie a bino mass m^ 143 GeV and corresponds to 
|2m^ — m^l < 3 GeV. For larger negative values of A^, rua increases, as can be seen from eq. (2.7), 
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Figure 1: Constraints in the mo, M1/2 plane for A = .1, Aq = —900 GeV, A^^ ~ —60 GeV and (a) tan/3 = 5, 
(b) tan/3 = 10. We show the region excluded by theoretical constraints or by LEP/Tevatron searches on 
sparticles (black), the region where a charged sfermion is the LSP (blue), the LEP limits from Higgs searches 
(red), the constraints from B-physics (pink) and from (g — 2)^ (violet) . The region allowed by the upper 
bound on the DM relic density fi/i^ is displayed in green. 

which means that rapid annihilation would be possible for a heavier bino LSP, that is a larger 
Ml/2- However, when is large it becomes increasingly difficult to rely exclusively on Higgs 
exchange to have efficient enough annihilation. Thus, for large negative values of Af^, the rapid 
annihilation region disappears. In addition, compatibility with 5a^ requires not too large values 
for mo, constraining further the allowed range of A,^ > — 100 GeV. On the other hand, if 

\Af^\ is too small (A^ ^ ~ 20 GeV), the a resonance band appears at small values of which 
are strongly constrained by chargino and Higgs searches at LEP at small tan/3, and by 6 ^ 57 as 
one increases tan/3. While determines the value of ^f which the a resonance appears, the 
value of Aq steps in the Higgs and sfermion masses: if |^o| is too small (^0 ^ ~ 500 GeV) the LEP 
Higgs bounds become too strong and the whole mo, M1/2 plane is excluded. Similarly, if \Aq\ is 
too large {Aq < — 1.5 TeV), the sfermions become too light, excluding the mo, M1/2 region where 
6a^ is within the 2a limits. Rapid annihilation of the bino LSP through the light scalar doublet 
h can also occur at low values of M1/2 ~ 130 GeV, as in the CMSSM. The low M1/2 region is 
however strongly constrained as explained above. Finally, as in the CMSSM, we find a bino-stau 
coannihilation band for values of mo just above the stau LSP forbidden region (at small mo) as 
well as a narrow bino-stop coannihilation region at small values of mo and M1/2 just above the stop 
LSP forbidden region (at M1/2 ~ 250 GeV). The upper limit on a^p from CDMS [40] is satisfied 
over the whole mo, M1/2 plane. 

At large values of tan/3 = 50, when ^Oi < 0, DM annihilation mechanisms are bino LSP 
annihilation through Higgs exchange, bino-sfermion coannihilation (at small mo) and higgsino LSP 
annihilation in the focus point region (large mo) . This is illustrated in fig. 2 for A^, = —60 GeV 
and (a) Aq = —900 GeV, (b) Aq = —1.5 TeV. The possible Higgs resonances are again the light 
scalar doublet h at low M1/2 ~ 130 GeV, just above the chargino exclusion limit from LEP or 
the pseudoscalar singlet a for slightly larger values of M1/2 ~ 200 GeV. The B-physics observables 
(predominantly b — > 57) restricts the small mo, M1/2 region, whereas the {g — 2)^ constraint 
restricts the large mo, M1/2 region eliminating most the focus point region. What is left of the 
focus point region by the {g — 2)^ constraint at large mo is excluded by B-physics constraints (here, 
predominantly Bg — > ^jl^ ji^ and — > T~^Vr). The LEP constraints on the Higgs sector excludes a 
large region of the parameter space, up to mo ~ 1.7 TeV for Aq = —900 GeV, eliminating all the 
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Figure 2: Constraints in the mo, plane for A = .1, tan/3 = 50, = —60 GeV and (a) Aq = —900 GeV, 
(b) Aq = —1.5 TeV. Same color code as in fig. 1 with in addition the limit from CDMS on a^p for {a-rrN, co) = 
(55,35) MeV (black dots) and (45,40) MeV (black dashes). 

bino-stau coannihilation region and leaving only a small region of resonance annihilation. However, 
over the whole region excluded by Higgs searches in fig. 2(a), the lightest Higgs doublet mass is just 
below the LEP limit (114.5 GeV). Hence, this constraint is easily relaxed by taking a larger value 
of l^ol) (^0 = —1-5 TeV in fig. 2(b)) or by assuming a 3 GeV theoretical uncertainty on the Higgs 
mass. Increasing |^o| also somewhat strengthtens the {g — 2)^ constraint by increasing the mixing 
in the smuon sector. The influence of Aq is most noticeable in the focus point region when is not 
too large since the smuon mixing is proportionnal to A^ — ^ tan/3. In the region where the heavy 
Higgs doublet H is light, the predictions for DM direct detection rate can be large. The CDMS 
limit thus constrains small values of Some area of the parameter space is excluded even when 

being conservative and fixing {(Tt^NjC^o) = (45,40) MeV, a choice that gives a small coefficient for 
the s quark content in the nucleon. 

For large tan/3 = 50 and Aq, A^ > 0, all Higgs states, except the scalar singlet s, can contribute 
significantly to the bino LSP annihilation. In fig. 3, we take Aq = 900 GeV and (a) A^ = 60 GeV, 
(b) Ai^ = 240 GeV. The possible Higgs resonances are: (i) the light scalar doublet h just above 
the chargino exclusion limit from LEP at M1/2 ~ 130 GeV; (ii) the pseudoscalar singlet a for 
Ml/2 ~ 200 GeV when A« = 60 GeV and M1/2 ~ 700 GeV when A^ = 240 GeV; (iii) the heavy 
scalar /pseudoscalar doublets H/A in the wide diagonal band. In addition, for such large values of 
tan/3, as in the CMSSM, one finds a region at large mo where the LSP has a significant higgsino 
component and therefore annihilates efficiently. Finally, at small tjiq, one still has a thin band 
where the bino LSP coannihilates with the stau LSP along the excluded region where the stau 
is the LSP. Note that regions at small niQ, M1/2 are in conflict with the LEP Higgs constraints 
as well as with the Bg — > /x^^^ and/or — > r'^i'r constraints (recall that for Aq > the 
b ^ s'j constraint is less effective). Regions at large niQ and small M1/2 are also excluded by 
the Bg fJ'^fJ'~ and/or i?+ T~^Ur constraints since the heavy neutral/charged Higgs doublets 
are rather light. Note that as one increases mo for a fixed value of M112, say 200 GeV, the 
charged Higgs mass mij± ~ rriA decreases until Br{B~^ — > t~^v-j-) drops below the allowed range 
for tuq ~ 2 TeV because of the destructive interference of the contribution with the standard 
model contribution. When niQ is increased further mu± continues to decrease such that 
Br{B^ — > T~^iyr) is large and dominated by the contribution. Hence, there is a narrow region 
where this constraint is satisfied, although for that region the branching ratio for Br{Bs — > 
is too large (in fig. 3 only the combined exclusion from B-physics observables is displayed). Then, 
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Figure 3: Constraints in the mo,Mi/2 plane for A = .1, tan/3 = 50, Aq — 900 GeV and (a) — 60 GeV, 
(b) Ai^ = 240 GeV. Same color code as in fig. 1 and fig. 2. 

as one increases rriQ, mfj± ~ uia decreases until it falls below the LEP limit for mo ^2.4 TeV and 
eventually becomes even negative for larger values of niQ. The constraint on {g — 2)^ excludes two 
distinct regions: at small mo, M1/2 the prediction for <5a^ is too large, and at large mo, M1/2 it 
is below the preferred range. The direct detection rate exceeds the limit from CDMS both when 
the LSP has a large higgsino component (the large mo focus point region) and/or when the heavy 
Higgs doublet H is light (at small Mi/2)- 

Setting the top quark mass at the minimum of the 2a range, nit = 171.2 GeV, one finds that 
a much larger area of the mo, plane satisfies the cosmological upper bound on the relic 

density when ^Oj^k > 0. This is because decreasing mt, all other parameters fixed, one obtains 
both a lower value for /x, hence a larger higgsino fraction for the LSP, as well as a lighter heavy 
Higgs doublet sector. The focus point and the H/A funnel regions of fig. 3 then merge in a large 
allowed band. This effect was already observed in [30], fig. 3(b) for a different choice of parameters. 
Furthermore, since is smaller, the smuon mixing is also smaller and the {g~2)^ constraint is less 
severe, especially in the focus point region. Even though the B-physics and LEP Higgs constraints 
are slightly more severe and a LSP with a larger higgsino fraction necessarily implies a larger direct 
detection rate, often exceeding the CDMS limit, one finds a larger area of the mo, plane that 
satisfies all constraints. Conversely, increasing mt = 174 GeV implies that the focus point region 
is pushed at higher values of mo and the H/A funnel appears at higher values of separating 
further the two DM allowed bands of fig. 3 . Similarly, B-physics and LEP Higgs constraints are 
slightly less severe when rat is larger, while {g—2)^ constraints are strengthened, moving the allowed 
area between these constraints in fig. 3 to smaller values of mo, -^1/2- The only remaining regions 
passing all constraints in this case are the bino-stau coannihilation region as well as the three Higgs 
funnels - the heavy doublet H/A, the light doublet h and the pseudoscalar singlet a, the latter two 
being however disfavoured by the CDMS constraints. Varying mt when Ao,Ak < has a similar 
effect on the location of the focus point region and on the phenomeno logical constraints. Yet in 
this case, there is no H/A funnel and the region in the mo, M1/2 plane that satisfies all constraints 
is only shifted, but not much larger for mt = 171.2 GeV. 

4.2 Small A: singlino LSP 

The singlino LSP can be found only for A ^ 1, this means A = .01 in our selected scans (qualita- 
tively, the results are similar for smaller values of A). In this case, we have an effective MSSM with 
an almost decoupled singlet sector and we do not expect to have singlet Higgs resonances. Thus the 
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Figure 4: Constraints in the toq, M1/2 plane for A = .01, tan/3 = 10 and (a) Aq = —30 GeV, A^, ~ —50 GeV, 
(b) Aq = 300 GeV, A^, = 50 GeV. Same color code as fig. 1, with in addition: the singlino LSP region with 
relic density above the WMAP limit (grey) and below WMAP limit (cyan). 

only mechanism that can provide the correct relic density for a singlino LSP is coannihilation, or 
more precisely self-annihilation of the NLSP. This works most efficiently when the LSP and NLSP 
are well below the TeV scale. We therefore expect to find WMAP allowed regions for choices of 
input parameters that predict a singlino LSP at small values of mo where the stau is light, or at 
large mo and tan/3, in the focus point region, where the higgsino is light. In the small singlet mixing 
limit (A <C 1) the singlino mass is given by mg = 2i/ which, according to eq. (2.3), implies 

m^ = 2{B-Ax). (4.1) 

Hence, the singlino mass depends mainly on mo, M1/2 (through B, see eq. (2.4)) and Aq (through 
A\ computed by integration of the RGEs) but is mostly independent of A^. It then suffices to take 
= iSO GeV (depending on the sign of Aq) in order to have positive mass squared for the Higgs 
singlet states. 

For intermediate values of tan/3 = 5 or 10, there is no focus point region and we can restrict 
ourselves to small values of mo where the singlino-stau coannihilation is expected to work out. In 
this region of parameter space however, it is difficult to satisfy both the LEP Higgs constraints, 
which exclude small values of the {g — 2)^ constraint, which exclude large values of 

Mx/2- For tan/3 = 5, we have checked, by fixing mo = and scanning over the only two remaining 
parameters, M1/2 and Aq, that the excluded regions overlap, leaving no possibility. We have 
repeated the same exercise for tan/3 = 10 and found that for Aq = —30 GeV, = —50 GeV 
and Aq = 300 GeV, = 50 GeV there are singlino LSP regions with the correct amount of DM, 
satisfying both the LEP Higgs and the {g — 2)^ constraints, see fig. 4. Note that the precise value 
of Aq is crucial for these plots. 

For Ak < (fig. 4(a)), the singlino mass, eq. (4.1), is positive (see eq. (2.8)) and grows with mo 
and Ml/2- Therefore, one finds a singlino LSP region at small Mi/21 below the stau LSP forbidden 
region. At the frontier of the two regions, one finds a band where the singlino LSP coannihilates 
with the stau NLSP, giving the correct amount of DM. For smaller values of Mi/2: the singlino LSP 
coannihilates with the stop NLSP. Large values of M1/2 ~ 650 GeV are excluded by the {g — 2)^ 
constraints, while the LEP Higgs constraints exclude small values of ^ 500 GeV. In the small 
region that passes the latter constraints near Mi 12 = 300 GeV, the light doublet h and the singlet 
s are almost degenerate in mass, and mixing effects push rrih just above the LEP limit. The LEP 
Higgs constraints can however be relaxed if one chooses a smaller value of A (e.g. A = .001) or a 
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Figure 5: Constraints in the mo, A/1/2 plane for A = .01, tan/3 = 50 and (a) ^0 = —600 GeV, = —50 GeV, 
(b) Ao = GeV, = 50 GeV. Same color code as fig. 4. 

larger value of mt, both solutions having the same effect of increasing ruh. Alternatively, one could 
allow for a 3 GeV theoretical uncertainty on m^, in which case the LEP Higgs excluded region 
reduces drastically. 

For > (fig. 4(b)), the singlino mass is negative and grows with rriQ and Mi^2- Therefore, 
the singlino LSP is found at large values of below the excluded region where the pseudoscalar 

singlet mass squared becomes negative (occurring at ~ 1-2 TeV). At the frontier with the 

forbidden stau LSP region, one finds a region where the relic density of the singlino LSP is below 
the WMAP limit. For our choice of ^0 = 300 GeV, this region satisfies in addition both the 
LEP Higgs and the {g — 2)^ constraints. However, the same remarks concerning the LEP Higgs 
constraints as in the case A^^ < apply here. In all cases, the mass difference between the singlino 
and the stau must be < 3 GeV for the coannihilation mechanism to be efficient enough. 

For large tan/? = 50, i? ~ (see eq.(2.4)) and the singlino mass is ~ — 2^a depending mainly 
on Hence, for large values of the bino is heavy and the singlino is the LSP. In addition, 
both the LEP Higgs and the (^f — 2)^ constraints are much weaker for large values of tan/?. In fig. 5 
we present the results for ^0 = —600 GeV, A,^ = —50 GeV and Aq = 0, A^ = 50 GeV. In both cases, 
the singlino is the LSP for M1/2 ^ 700 GeV and its relic density is below the WMAP upper bound 
when it coannihilates with the stau NLSP at low mo. For M1/2 ~ 700 GeV, as in the CMSSM 
the bino LSP also coannihilates with the stau NLSP at low mo or annihilates through the light 
Higgs doublet h resonance when M1/2 ~ 130 GeV. The latter region is however mostly excluded by 
the B-physics and LEP Higgs constraints. As in the CMSSM, we find a WMAP allowed band at 
large mo, where the LSP has a significant higgsino component and therefore annihilates efficiently. 
Most of this band is however excluded by CDMS, even when taking the conservative limit. At 
the frontier between this region and the singlino LSP region, we find a thin WMAP allowed band 
where the singlino LSP coannihilates with the higgsino. Yet, this region is incompatible with the 
(5 ~ 2)^ constraint. For Af^ > 0, the heavy Higgs doublet H/A resonances can also contribute to 
the annihilation of the bino LSP in a large diagonal band, as in fig. 3. However, in the singlino LSP 
region, when M1/2 ^ 700 GeV, this mechanism is not sufficient to bring the relic density below the 
WMAP upper bound, except in a thin band at the frontier between bino and singlino LSP regions 
where the singlino LSP coannihilates with the bino NLSP which in turn annihilates through the 
H/A resonance. 
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Figure 6: Contours of a^p (in pb) in the mo, plane for A = .1,Aq = —900 GeV, = —60 GeV and 

(a) tan/3 = 5, (b) tan/? = 10. 



5 Prospects for dark matter direct detection 

In sec. 4 we have seen that direct detection experiments (CDMS) already constrain the parameter 
space of the CNMSSM at large values of tan/3, especially in the focus point region at large tjiq. 
Here we present in more details the predictions for the SI cross sections for elastic scattering 
of neutralinos on protons in scenarios with an extra Higgs resonance. Cross sections for direct 
detection of a singlino LSP are always far too small to be within the present or future experimental 
reach. Recall that experimental sensitivities a^j, ~ 10~^ — 10~^ pb are expected within a few years 
and that the planned ton-scale detectors such as Warp, Xenon, Eureca [95] or SuperCDMS [96] 
will reach a sensitivity around 10"^*^ pb within a decade. We will see that for large regions of 
the parameter space of the CNMSSM the predictions for SI direct detection cross sections are well 
within reach of the ton-scale detectors. In many cases signals are to be expected much before that. 
We take here the micrOMEGAs default values, (cr7rAr,o"o) = (55,35) MeV to determine the quark 
coefficients in the nucleon. Almost one order of magnitude variation in the prediction can result 
from diffrent choices for the quark coeficients [87, 97, 98]. 

At small to intermediate values of tan/3, a^j, is dominated by the light Higgs doublet h exchange 
although both the heavy Higgs doublet H and squarks exchange diagrams also contribute. Contour 
plots of a^p for the cases considered in the previous section, X = .1, Aq = —900 GeV, Af^ = —60 GeV 
with tan/3 = 5 and 10, are displayed in fig. 6. (In the regions where no contours are displayed, 
either the stop or the stau is the LSP). In the region favoured by all constraints (see fig. 1), the SI 
cross section is predicted to be between 5(2) and 10 x lO"^'^ pb for tan/3 = 5(10) with a LSP mass 
between 70 and 200(250) GeV (here the bino LSP mass is simply given by ~ AlMi/2)i within 
reach of the ton-scale detectors. Note that typically a^j, decreases with M1/2, as both the mass of 
the heavy Higgs doublet and of the squarks increase with M|/2- In the bino-stau coannihilation 
region at small mo, large M1/2 (which is however disfavoured by the {g — 2)^ constraint) the SI 
cross section drops below the sensitivity of ton-scale detectors as in the CMSSM. 

For tan/3 = 50, a^p is dominated by the heavy Higgs doublet H exchange since the couplings 
of H to quarks and to the LSP are tan/3 enhanced. This enhancement compensates for the mass 
suppression factor mf^/rnjj. Cross sections near or above the present limit are obtained for small 
Ml/2, see fig. 7(a) for A = .1, Aq = —900 GeV and Af^ = —60 GeV. In the allowed region where 
the LSP annihilates through the pseudoscalar singlet a resonance (M1/2 ~ 200 GeV, see fig. 2) we 
find a^p ~ 2 x 10"*^ pb. This region should soon be probed experimentally. A large negative value 
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for Aq = —1.5 TeV induces larger heavy Higgs masses, thus a smaller SI cross section, see fig. 7(b). 
In this case, the pseudoscalar singlet funnel corresponds to a^p ^ 2.5 x 10~^ pb. In both plots of 
fig. 7, regions where no contours are displayed correspond either to regions where no solution to 
minimization eqs. (2.4) is found (large mo), the LSP is a sfermion (small mo). In the bino-stau 
coannihilation allowed region (small mo, 600 GeV < Mx/2 ^ 1 TeV, see fig. 2) the SI cross section 
is found to be barely within the reach of the ton-scale detectors in both cases. On the contrary, as 
mentioned in sec. 4, the predictions for a^p are already above the present experimental bounds for 
most of the focus point region (large mo, small M1/2) unless one allows for reduced coefficient of 
the s quarks in the nucleon. 

Larger cross sections can be found in the bino-stau coannihilation region for different choices 
of ^0) since the mass of the heavy Higgs doublet H as well as its coupling to the LSP depend on 
Aq. Contours for the same parameter choice as in sec. 4, A = .1, tan/3 = 50, = 900 GeV and 
Ak = 60 GeV, are displayed in fig. 8. In the bino-stau coannihilation allowed region (small mo, 
700 GeV < Ml/2 < 1.2 TeV, see fig. 3) a^p is now in the range 0.6—6 x 10"^ pb, i.e. within the future 




in„ [GeV] 

Figure 8: Contours of a^p (in pb) in the toq, M1/2 plane for A = .1, tan/3 — 50, Aq = 900 GeV and 
Af^ = 60 GeV (contours for A^ = 240 GeV are identical). 
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experimental reach. The pseudoscalar singlet funnel {M1/2 ^ 200 GeV, see fig. 3) now corresponds 
to a^p > 2 X 10~^ pb which is already excluded. The value of A,^ has no influence on the SI cross 
section since the dominant process is the heavy Higgs doublet H exchange and influences only 
the singlet Higgs a/s masses. Contours for Af^ = 240 GeV are therefore basically identical to those 
for A^ = 60 GeV. In this case however, the pseudoscalar singlet funnel {M1/2 ~ 700 GeV, see fig. 3) 
corresponds to ^ 3 x 10~^ pb which should be probed by the ton-scale detectors. Finally, as for 
the Ak <0 scenarios, most of the focus point region is already excluded by DM direct detection. 



6 Higgs searches at the LHC 

As mentioned in the introduction, the discovery of three neutral Higgs states with large mass 
sphttings would provide a strong evidence in favour of the NMSSM. Furthermore a determination 
of the masses of the neutral Higgs states, in addition to the measurement of the LSP mass and 
couplings, is important for a precise theoretical prediction of the DM relic density and/or direct 
detection cross section. Here we do not address the issue of how precisely masses of sparticles and 
Higgs bosons can be determined experimentally, however as a first step we explore the discovery 
potential for the neutral Higgs states of the NMSSM at the LHC. We concentrate on the regions 
in the parameter space allowed by all phenomenological constraints, which are: the h, a and H/A 
funnels and the stau coannihilation region (the focus point region is disfavoured by the (5 — 2)^ 
and DM direct detection constraints). 

The most relevant modes for detecting the NMSSM neutral Higgs bosons at the LHC are those 
that have been studied for the SM and the MSSM. These are: 

1) associated bbcj) production with (j) — >■ t~^t~ 

2) M ^ ^ ^ 77 

3) M ^ ^ ZZ(*^ 4 lept ons 

4) M ^ ^ VI7iy(*) ^ i+i-i^i? 

5) WW ^ ^ r+r" 

6) WW -^<P^ WW 

7) WW ^ ^ 77 

where £ = e, jj, and ^ stands for any neutral Higgs scalar or pseudoscalar state, the latter only for 

channels 1) and 2) though, as a pseudoscalar state does not couple to WW or ZZ at tree level. 
To compute the statistical significances {No- = S/ of all Higgs states in each channel, wc have 
followed the same procedure as in ref. [18]. For channels 2-7) this consists simply in multiplying 
the expected N^r in the SM by the reduced coupling of the considered NMSSM Higgs state to 
Icptons or gauge bosons as compared to those of the SM Higgs with the same mass. For channel 
1) the issue is more complicated as the experimental reach is usually given as 5a contours in the 
niA, tan/3 plane. In the MSSM, these contours correspond to associated production of the heavy 
doublet states H/A in equal proportion as: (i) ruA ~ ruH within the r+r" mass resolution {ii) 
BR{A T~^T~) ~ BR{H T^T^) ~ 10% in the considered mass range and {in) the bbA and 
bbH couplings arc equal and scale as tan/3. The last point explains why this channel is especially 
important in the MSSM at high tan/? but not considered for the SM Higgs boson. In addition, 
experimental cuts designed to single out the associated production (mainly 6-tagging) allow us to 
neglect the contribution of gg H/A to the signal. As a result, the net signal rate along the 5(T 
contours is twice that for bbH or bbA alone. Thus, N^j =2.5 would be achieved for bbH or bbA 
along these contours were niA and mn widely separated. Defining the excluded value of tan/3 = t 
as a function of ruA along the contour plots, we compute 
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where (/) stands for any NMSSM neutral Higgs scalar or pseudoscalar state and 6^ is its reduced 
coupling to b quarks as compared to that of the SM Higgs with the same mass. Finally, in each 
channel we have combined the signals for degenerate Higgs states with a mass resolution of 15% 
for channel 1), 1% for channels 2-4) and 10% for channels 5-7). In ref. [18] high luminosity {C = 
300 fb""*^) optimistic results from CMS and ATLAS were used. Here, we take the updated and more 
realistic low luminosity {C = 30 fb~^) predictions from CMS, summarized in fig. 11.35 for channel 
1) and fig. 10.39 for channels 2-7) of the CMS-TDR [99]. 

In order to explore the LHC reach for NMSSM Higgs searches, we have plotted contours of the 
maximal significance for each Higgs state in the plane mo, M1/2 with the same parameter choices as 
in sec. 4-5. Rather than showing a large number of plots, we will now discuss some general properties 
and eventually present a few selected case studies. First, let us remark that for all parameter choices, 
the maximal significance for the lightest Higgs doublet h is always Nu{h) ~ 4, but it hardly reaches 5 
in a single channel. This situation (similar to that of the MSSM) means that a high luminosity run, 
improved detection techniques or channel combinations will be necessary in order to exhibit a clean 
5(7 signal for the light h at the LHC. The best channels for h discovery are either 2) gg ^ h ^ 
or 5) WW — >■ /i — > t'^t~. However, other channels as 3) gg ^ h ^ ZZ^*^ — >■ 4 leptons and 
7) WW ^ /i ^ 77 can yield a significance ~ 1 — 3 (we have not combined significances in different 
channels for each Higgs state). The heavy doublet states H/A are always nearly degenerate and 
the LHC will not be able to see them separately. The best channel for these states is 1) associated 
bbH/A production with H/A t~^t~ . This channel covers masses up to mujA ~ 800 GeV and is 
enhanced by a factor tan^/3. Hence, a ba signal shall be seen at the LHC provided that Mi 12 is 
not too large and tan/3 is large. In our selected scans, we indeed found that for tan/3 = 5 or 10 the 
significances for H/A are always small and N„{H) = N^{A) = 5 is reached only for tan/3 = 50 and 
Ml/2 iiot too large. 

In the singlet sector, we found that the scalar s is usually too heavy and/or too decoupled to be 
seen at the LHC. This is because if the singlet s is light and A is not too small {ie the doublet /singlet 
mixings are substantial) then mixing effects between h and s always lead to a light state in the 
Higgs spectrum, already excluded by LEP constraints (note that this is not the case for a light 
pseudoscalar singlet a which does not mix with the light h). On the other hand, if the singlet s 
is light and A is small (ie the doublet/singlet mixings are negligible) then one obtains an effective 
MSSM with a decoupled singlet sector. One can still have differences between such a scenario and 
the true MSSM if the singlino is the LSP (cf. sec. 4.2). However, the couplings of the singlet Higgs 
states to matter fermions and gauge bosons are then so small that it would be impossible to see 
them at the LHC (we will not discuss here the possibility of observing a singlino LSP at the LHC). 
The significances for the scalar singlet s are therefore always small. The only exception is when 
it is nearly degenerate with the heavy doublet states H/A. In this case, the significances for the 
three states H, A and s are equal and can be large, although the net contribution of the singlet s 
to the total significance is always small and the LHC would see only one peak. 

The pseudoscalar singlet a can play an important role in DM reHc density computation, as 
we have seen in sec. 4.1. Contrary to the scalar singlet s it can be relatively light without being 
excluded by LEP. If it is not too decoupled (ie if A is not too small) one can then expect to be 
able to see it at the LHC. For the selected scans of sec. 4-5, we have checked that for A = .01 
the significances for the singlet a are always small. For A = .1 we have found some cases where 
these significances are large enough to be able sec something at the LHC. These cases correspond 
to tan/3 = 50 only as the main discovery channel for the singlet a is 1) associated bba production 
with a T^T~ , which is tan/3 enhanced. We will now study in details a few cases, summarized in 
table 1. 

Let us start with ^0 = —900 GcV, A,^ = —60 GeV as in fig. 2(a). After imposing all phenomeno- 
logical constraints, the only region allowed is the a funnel for ^ 2 TeV and M1/2 ^ 200 GeV. 
Case (1) in table 1 corresponds to a point in this region. The significance for the pseudoscalar sin- 
glet in this case is very small (iVo-(a) = .1). Here, the mass of the scalar singlet is = 567 GeV, 
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Case 


^0 




mo 


Ml/2 






ruh 




rriH 


rriA 


N^{H/A) 


•ma 




(1) 


-900 


-60 


1900 


207 


.129 


87 


114 


4.1 (2) 


528 


527 


13 (1) 


183 


0.1 (1) 


(2) 


-1500 


-60 


1600 


130 


.133 


56 


115 


4.2 (2) 


434 


441 


17(1) 


191 


1.0 (1) 


(3) 


900 


60 


1625 


155 


.130 


61 


114 


4.4 (5) 


296 


309 


26 (1) 


190 


3.9 (1) 


(4) 


900 


60 


1630 


305 


.035 


123 


114 


4.4 (5) 


315 


298 


27(1) 


200 


5.9 (1) 


(5) 


900 


60 


900 


650 


.133 


273 


114 


4.2 (5) 


476 


482 


19 (1) 


251 


0.8 (1) 


(6) 


900 


60 


475 


740 


.110 


311 


114 


4.2 (5) 


420 


430 


23 (1) 


253 


1.7 (1) 



Table 1: Higgs masses and statistical significances at the LHC with C = 30 fb for selected points. 
In all cases tan/3 = 50, A = .1 and mt = 172.6 GeV (except case 4 where mt = 171.2 GeV). The 
channel where the significance is maximum is indicated in parenthesis. Dimensionful parameters 
are in GeV. 

close enough to niH, tua to be within the t'^t~ resolution. Hence N(j{s) = N„{H) = N^j^A) = 13, 
although the LHC would see only one peak. The significance for the light doublet N(j{gg ^ h ^ 77) 
is below 5. However, Na{gg ^ ZZ^*'^ 4 leptons) = 2.7, Na{WW h ^ r+r") = 4.0 and 
N^{WW h ^ 77) = 2.2, so that combining the signals from different channels would yield a 
global significance above 5. The situation for the light h is similar in all our case studies. 

For Aq = —1500 GeV, Af^ = —60 GeV as in fig. 2(b), there are 3 regions allowed by all 
phenomenological constraints. In the a funnel, as well as in the stau coannihilation region, one 
finds very small significances for the singlet a, Ncr{a) ~ .1. However, in the h funnel one finds slightly 
larger significances. This is illustrated by case (2) in table 1 where N^{a) = 1. Higher luminosity 
or improved detection techniques could eventually lead to a 5cr signal for the pseudoscalar singlet 
in this scenario. 

Next, we consider the case where Aq = 900 GeV, A^ = 60 GeV as in fig. 3(a). Here, the allowed 

regions are the h, a and H/A funnels, the stau coannihilation region and the focus point region 
(although the latter is usually disfavoured by the (g — 2)^ and DM direct detection constraints). 
Case (3) in table 1 corresponds to a point in the h funnel where all significances are quite large, 
making it possible to distinguish the NMSSM from the MSSM. Even larger significances for the 
singlet a (Ncj(a) > 5) are found at smaller values of M1/2, although these are disfavoured by the 
T^Vr constraints. Note that with a smaller value for the top quark mass the doublet /singlet 
{A/ a) mixing is increased and the significance for the singlet a is already above 5 as can be seen in 
case (4) of table 1. This point passes all experimental constraints. However, the bino annihilation 
cross section is large due to the singlet a resonance and the DM relic density is below the range 
allowed by cosmological data. In such a scenario, one would then have to rely on a different 
mechanism to explain the abundance of DM in the universe. Finally, case (5) sits in the H/A 
funnel and case (6) in the stau coannihilation band. In both cases, we find substantial significances 
for the singlet a, although always below 5a. As for case (2), a higher luminosity or improved 
detection techniques could eventually lead to a 5a signal for the singlet a in these scenarios. Larger 
values of A^, eg A^ = 240 GeV as in fig. 3(a), yield a heavier singlet a and therefore smaller 
significances. 

7 Conclusion 

We have reexamined the DM favoured scenarios in the CNMSSM taking into account constraints 
from B-physics observables and from the (5— 2)^. Although the allowed parameter space is strongly 
constrained by these observables we confirm our previous results : scenarios where the LSP annihi- 
lates near a pseudoscalar Higgs singlet resonance or where the LSP is mainly singlino are allowed 
in addition to the MSSM-like scenarios such as bino-sfermion coannihilation, bino annihilation 
through a Higgs doublet resonance or mixed higgsino LSP. 

B-physics observables constrain mainly the scenarios at large tan/3 with b ^ sj being more 
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effective in tlie low mo, M1/2 region when Aq is negative while fi~^fj~ and t^Uj- being 

more powerful constraints when ^0 > and/or at large mo- The most powerful constraints on the 
parameter space are however obtained from the {g — 2)^. Indeed, explaining the deviation from 
the SM by SUSY contributions requires light sfermions and charginos, especially at low values of 
tan/3. The focus point region at large mo as well as the bino-sfermion coannihilation region at large 
Ml/2 s-re therefore strongly constrained. This feature is not specific to the CNMSSM and is also 
observed in the CMSSM. A significant portion of the Higgs (doublet or singlet) funnel regions is 
also constrained by [g — 2)^. 

The recent upper limits on the neutralino proton elastic scattering cross section constrain mainly 
the focus point region at large tan/3. When the LSP is partially higgsino and tan/3 is large, the 
predictions from rise very fast with the decrease of the mass of the heavy Higgs doublet so 
one will easily probe all the allowed parameter space in the near future. This is so even when 
taking into account the large uncertainties in the theoretical predictions. In fact we have shown 
that future direct detection experiments have good prospects for probing the parameter space of 
the CNMSSM with a dominantly bino or higgsino LSP. On the other hand, at small values of 
tan/3, or when the Higgs and squarks are in the TeV range, one will have to wait for the large 
scale detectors. This statement however is strongly correlated with the lower limit on (5a^ we 
have imposed which disfavours the scenarios at large values of M1/2 with little prospects for direct 
detection. The large theoretical uncertainties in the computation of the direct detection rate is an 
issue. In particular a better determination of the quark coefficient in the nucleon could significantly 
reduce the uncertainty. Finally the characteristic NMSSM scenario with a singlino LSP is far beyond 
the reach of large scale detectors. Similar conclusions were reached in an analysis performed in the 
context of another constrained version of the NMSSM where A^ = Aq [100] . 

In the near future with the onset of the LHC, improved sensitivity on B-physics observables and 
improved direct detection results, prospects for discovery of supersymmetry are good. If a clean 
signal of physics beyond the standard model is observed, the issue will be to establish whether or 
not the CNMSSM (or even only the NMSSM) is the correct scenario. What distinguishes the model 
from the MSSM is the more elaborate neutralino and Higgs sectors. We have analysed the potential 
of the LHC to unambiguously identify the NMSSM by observing at least 3 neutral Higgs states. We 
have found that this is possible only for large values of tan/3 ~ 50 and A ~ .1 where one can benefit 
both from the enhanced couplings of HjA to hh and r+r" and from substantial doublet /singlet 
mixings. The best signals are found when the heavy doublets H/A and the pseudoscalar singlet a 
are at most a few hundred GeV's (that is at low M1/2 or mo). This region is also the one that will 
be efficiently probed by direct detection searches. At lower values of tan/3, there is little prospects 
of discovery at the LHC for the Higgs singlet states. This means that in the scenarios where 
the pseudoscalar singlet a contributes significantly to the bino annihilation, the only indication in 
favour of the NMSSM would be that the MSSM prediction for the DM relic density is much above 
the cosmological measurements. 

Note also that a very Hght pseudoscalar singlet a can appear in the NMSSM and it is not 
completely excluded by LEP constraints nor by recent T radiative decay constraints [21-24]. It 
could even explain the small 2.3o" excess in the e+e^ Zbb channel at LEP if the light doublet h 
decays mainly as h ^ aa [21, 101-103]. In this case, new detection channels might reveal necessary 
in order to reinstate the 'no-lose' theorem for NMSSM Higgs discovery at the LHC [18, 19,46-51]. 
Finally, the observation of a singlino at the LHC or at the ILC would of course be another method 
to provide evidence for the NMSSM [52-57]. 
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Note added 

A recent estimate of the hadronic contribution to (g — 2)^ based on the latest data from BaBar 
indicates that 5a ^, the deviation between the measured value and the theoretical SM prediction, 
may be smaller than considered here [104]. If this result is confirmed, the Ja^u constraint in the 
large M1/2 region would disappear. 
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